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ABSTRACT 

Turbulent convection models are thought to be good tools to deal with the convec- 
tive overshooting in the stellar interior. However, they are too complex to be applied 
in calculations of stellar structure and evolution. In order to understand the physical 
processes of the convective overshooting and to simplify the application of turbulent 
convection models, a semi-analytic solution is necessary. We obtain the approximate 
solution and asymptotic solution of the turbulent convection model in the overshoot- 
ing region, and find some important properties of the convective overshooting: I. The 
overshooting region can be partitioned into three parts: a thin region just outside the 
convective boundary with high efficiency of turbulent heat transfer, a power law dis- 
sipation region of turbulent kinetic energy in the middle, and a thermal dissipation 
area with rapidly decreasing turbulent kinetic energy. The decaying indices of the tur- 
bulent correlations k, u'^T' , and T'T' are only determined by the parameters of the 
TCM, and there is an equilibrium value of the anisotropic degree uj. n. The overshoot- 
ing length of the turbulent heat flux u'^T' is about lHk(Hk = IJ^D- m. The value 
of the turbulent kinetic energy at the convective boundary kc can be estimated by a 
method called the maximum of diffusion. Turbulent correlations in the overshooting 
region can be estimated by using kc and exponentially decreasing functions with the 
decaying indices. 

Subject headings: convection — diffusion — turbulence 
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1. Introduction 

Convective overshooting is an important physical process in the stellar structure and 
evolution. Phenomenologically, the acceleration of a fluid element is zero at the convective 
boundary, but its speed is not zero. It is able to go across the convective boundary into 
the dynamically stable zone. This phenomenon is called the convective overshooting. The 
convective overshooting transports heat and matter, and affects the structure and evolution of 
stars. A phenomenological theory of the overshooting was developed by Zahn(1991), which 
predicts an adiabatic overshooting region. However, Xiong & Deng(2001) pointed out that 
the turbulent velocity and the temperature are strongly correlated in Zahn's theory. Recently, 
Christensen-Dalsgaard et al.(2011) found that the convective overshooting only described by the 
turbulent convection models could be in agreement with the helioseismic data. 



The turbulent convection models (TCMs) are based on fully hy drodynamic r nome nt 



equations, and applied on investigating the convective overshooting(] 
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). There are two main difficulties restricting the applications of the TCMs. One is to solve 



the equations of the TCMs, which are highly non-linear and unstable in numerical calculations. 



The other is to incorporate the TCMs into a stellar evolution code. In general, solving the TCMs 
needs the parameters of the stellar structure(e.g. temperature T, density p, pressure P, radius 
r, luminosity L, and elements abundance vector), and solving the equations of stellar structure 
requires the temperature gradient V which is determined by the TCMs. Thus, in order to apply 
the TCMs, one must solve both the TCMs and the equations of stellar structure, which shows 
enormous difficulty. Although developing numerical technique is very important, getting an 
approximate solution of the TCMs is more interesting because an approximate solution helps to 
understand the physical processes and may significantly simplify the application of the TCMs. 
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Xiong(1989) found the asymptotic solution of his TCM in the overshooting region, the turbulent 
correlations being exponentially decreasing in the overshooting region. However, his solution of 



the heat flux u'^T' is not suitable near the convective boundary, and the initial turbulent kinetic 
energy is unknown so that the value of the turbulent correlations in the overshooting region 
actually can not be determined without numerical calculations. 



In this p aper, we investiga te the properties of the convective overshooting 



Yang's TCMdLi & Yang 



Yang & Li 



20071), which was tested in the solar convection zone (|Li & Yang 



by analyzing Li & 



2007 



2007|). We try to get a semi-analytical solution of the TCM in the overshooting region. 



We introduce the TCM in Section 2, investigate the properties of the overshooting in Section 3, 
and summarize the conclusions in Section 4. 



2. Turbulent Convection Model 



The closure assumptions of Li & Yang's TCM are (|Li & Yang 



2001, 



20071): the three-order 



moment terms are modeled with a gradient-type scheme; the dissipation rate e of the turbulent 



kinetic energy k is assumed to be local; the dissipation rates of the turbulent heat flux u[T' 



and the turbulent fluctuation of temperature T'T' are assumed to be determined by both the 
reciprocal timescale of the turbulent dissipation t^^ = | and the thermal dissipation one 
""^ — ^p-. According to those closure assumptions, fully hydrodynamic moment equations 



on th e quasi- steady ap proximation result in the complete equations of two-order moment 



terms dLi & Yang 



2007b: 



1 ^ ( C or^^W^'''''''^ 



1 d 



,k 



-.Y,{CsPr-^.^)=e + ^l^ 



dr 



T 



(1) 
(2) 



2 d / ^k—-du'^T 
^"/^^ — 

p^z (jf- y £ (j^ 



u',T' (3) 
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^ ^ ^Ceipr2-<<^') = -^(V - V..)<F+ 2Ce ( r + — ) TT (4) 



pr'^ dr \ e ^ ^ dr J Hp ^ \k pcp 

The temperature gradient is calculated as: 

V = Vn-^'-^<r (5) 



T he meaning of those equations and each term in them were described in previous 



works (|Li & Yang 



2007 



Zhang & Liii2009|) in detail. We simply introduce them here: 



Equations (1-4) describe the equilibrium(time-independent) structure of the radial kinetic 



energy u'^.u'^, the turbulent kinetic energy k, the turbulent heat flux u'^T' and the turbulent 



fluctuation of temperature T'T', respectively. On the left side of those equations, there is the 
non-local term(i.e. the diffusion term) of each turbulent correlation. On the right side, there are 
the local terms which describe the generation and the dissipation of each turbulent correlation. 

3 

In Eq.(l) and (2), e is the turbulent dissipation rate of k and e = ^ where / = aHp, and the 
second term on the right side is the generation rate of the kinetic energy due to the contribution of 
the buoyancy. The last term in Eq.(l) is the return to isotropy term which attempts to make the 
turbulent motion be isotropic. In Eq.(3), the first two terms on the right side is the generation rate 



of the turbulent heat flux u'^T', and the last one is the dissipation rate that comprises the turbulent 
dissipation and the thermal dissipation. In Eq.(4), the first term on the right side is the generation 



rate of the turbulent fluctuation of temperature T'T', and the last one is the dissipation rate. 
Meanings of other symbols are: Hp = — is the local pressure scale height, /3 = — (^^)p 

dlnT 
dlnP 



the expansion coefficient, Qr = — the radial component of gravity acceleration, V 



the temperature gradient in the stellar interior, Vad = {§^)s the adiabatic temperature gradient, 
A = the thermal conduction coefficient, cp = (|f-)p the specific heat, Ck the parameter of 
the return to isotropy term, (Cs, Ca, Cei) the diffusion parameters and (a, Ct, Ce) the dissipation 



parameters of turbulent variations(A;, u'^T', T'T'). 
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In Eqs.(l-4), overbars are only used in three turbulent correlations u'^u'^, u'^T' and T'T' . The 
other variations(density p and the temperature T, etc.) are all mean state quantities which should 
use overbars but we ignore them for convenience. 

Equation (5) describes the energy transport in the stellar interior by both turbulent motions(i.e. 
convection and overshooting) and radiation. Vr is the radiative temperature gradient. 



3. Theoretical analysis of TCM in the overshooting region 



In the previous work dZhang & Lil 120091) . we applied the TCM in the solar overshooting 



region and found some properties of the overshooting region: u'^T' < 0, V r < V < W ad, and 



the peak of T'T', which are similar to Xiong's(1985) and Xiong & Deng's(2001) works. In this 
section, we attempt to get semi-analytical solutions of the TCM. 

Some approximations are adopted to simplify Eqs.(l-5) in the overshooting region: 

Approximation I. Peclet number Pe > 1, where = That is | > which 

means the turbulent dissipation is much stronger than the thermal dissipation. This assumption is 
reasonable in most cases except for the region near the surface of a star or with very small k. 

Approximation 11. All variations, except the turbulent fluctuations, are thought to be constant 
because the turbulent fluctuations change much faster than others in the overshooting region. 

Approximation 111. Far away from the convective boundary, V ~ V r. This assumption is 



acceptable if the heat flux u'^T' is small. 
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3.1. Hirbulent heat transport in the overshooting region 



Defining the anisotropic degree uj — which is the ratio of radial kinetic energy to total 
kinetic energy, and applying Approximation II and Eq.(5), we can rewrite Eq.(3) to: 

/ Sti' T' \ T Ba \/k 

— UCnUjlVk-^j = -j^i^R - V„,)« + ^T'T' + [2uP, + ^(1 + P'^)] (6) 

In the last bracket in Eq.(6), Approximation I{Pf. >> 1) makes the dissipation term 
Ct(l + P-^)^<T' be ignorable. And, by using Eq.(5) and Approximation II, it is easy to find 
that the diffusion term is on the same order of the ignorable dissipation term: 



dr \ 



Therefore the diffusion term is also ignorable. Equation (3) is in local equilibrium: 

- ^ ( V - Vad)KK + ~ (8) 

lip 1 

In the overshooting region, the most important process is the diffusion of the kinetic energy. 



Thus, we ignore the diffusion of T'T'(i.e., setting Cei = 0). The solution of the TCM with 
Cei = can be thought as the zero-order solution of the TCM. 



Ignoring the diffusion of T'T' and the diffusion and dissipation terms of u'^T', using 
Approximations I & II, one can rewrite Eqs.(l-4) as: 



3 



a:) = - 1)(- 5)- + f - (9) 



2c4(..i|,4 + f<r (10) 



0=-f^(V-VJa;*:+^T'r' (11) 
tip 1 
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2T £ 

■—{V-Vad)KT' + 2Cej^T'T' 
rip K 



(12) 



Equation (9) results from Eq.(l) and (2), describing the equilibrium structure of the 
anisotropic degree u. The left side is the diffusion of u. The first term in the right side is the 
dissipation rate due to return to isotropy term in Eq.(l). The last term is the generation rate of to 
due to the buoyancy. 



Equations (11) and (12) show: 



T kl 



{V -Vad){KT' + 2CeiO — — 

I3g I 



The solution is u'^T' = —2Ce£UJj^ or V = Va^. The latter is equivalent to u[.T' 



(13) 



■^-^C^ad — Vj?). Because m'^T' is close to z ero near the convective boundary and gradually 



decre ases far away from the convective boundary (|Xiong 



1989; 



Xiong & Deng 



2001 



Zhang & Li 



20091), the physically acceptable result is: 

T A 



Max{-^ — {Vad - Vk), -2CeU^^} 
Hp pcp Pg I 



(14) 



Equation (14) shows that there is an adiabatic stratification zone in the overshooting region 
in the case of Cei = 0. In order to investigate the property of heat transport in the overshooting 
region, we must know the length of the adiabatic stratification zone. It is found in Eq.(14) that the 

3 

boundary of the adiabatic stratification is the location where -^-^C^ad — Vi?) = 2CeUJj^^. 
Solving the equation of uj is not easy because it is nonlinear. However, this problem is avoidable. 
Turbulent motions are isotropic when oj = |. In the convection zone, w > | because the buoyancy 
boosts radial turbulent motion. In most part of overshooting region, cu should be less than | 
because the buoyancy prevents radial turbulent motion. Therefore cu should be not far away from 
I near the convective boundary. Further more, taking a; as a constant, one can rewrite Eq.(lO) as: 



(15) 
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Substituting Eq.(14) into the above equation, one can get the approximate solution: 



3 i / 3 r — Tn 

« ^c(^^P{-\l (16) 



if J^^(Va.-V«)<2Ce(^j4,and: 



(17) 

In Eq.(16), point C, which is the convective boundary where Vad — Vr, is set to be the initial 
point, kc and being k and r here. The contribution of the buoyancy temi(i.e. the last term 
in Eq.(15)) is ignored in obtaining the solution Eq.(16). In the deep convection zone, turbulent 

3 

motions are almost in local equilibrium, thus the ratio of —^u'^T' to ^ is about 1. However, 
near the convective boundary, buoyancy is about zero, meanwhile the diffusion of k dominates. 
Those make the ratio be much less than 1. Therefore the buoyancy term is ignorable. 

In Eq.(17), point A, where k = kA and r = r^, is the boundary of the adiabatic overshooting 
region. In the region beyond point A, the ratio of —^u'j.T' to ^ is 2CeOJ which is on the order of 
1, thus the buoyancy term remains. 

The exponentially decreasing function of k is due to the fact that there is no generation in the 
overshooting region. Contrary to the situation in the convection zone, the buoyancy dissipates k 
because it prevents the radial motion of fluid elements in the overshooting region. The distribution 
of k results from the equilibrium between the diffusion and the dissipation, k should decrease 
faster if the buoyancy is as effective as the turbulent dissipation, which is found by comparing the 
exponential indices of Eq.(16) and (17). 

3 

The location of point A is determined by -^-^{^ad — Vr) — 2CeUJj^^. Using Eq.(16), 
we get a property of point A: 



12 ( / 3 .TA-rc,. 1 a/3gX 
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The relation between va and Vr^a is needed in order to solve this equation and to locate 
point A. Near the convective boundary, there is: 

\^ad - Vr,a\ ^ Vad\x{lnPA - lnPc)\ = Vadlxl^ (19) 

Hp 

where lad — \ fA — fc\'^^ the length of the adiabatic overshooting region. Pa and Pc the pressure 
at point A and C, and x = which is approximately a constant. 

Substituting Eq.(19) into Eq.(18), one finds: 



A , I 3 lad. 1 a^gX lad .^f.. 

k^exp[ — J = — VadXTT ^^^^ 

ay 'iCsOJ Hp 2CeW pcp Hp 

lad can be worked out if kc is known. In the deep adiabatic convection zone, turbulent 

3 

diffusion is ignorable, and the localized TCM shows /c^ocai ~ °'^9K^^r-^ ad) ^^gg Appendix A). 
However, kc can not be estimated as that because — Vad thus kiocai = at the convective 
boundary. Actually, the turbulent diffusion of k is effective near the convective boundary, and kc 
is determined by the diffusion. We can estimate kc by a simple approach which will be referred to 
as the maximum of diffusion hereafter. Setting point B at where the diffusion becomes dominative 
in the convection zone, we get the relation between kc and ks by solving Eq.(15): 



k-^ = klexp{-J-—-\'-^—^\) (21) 



ACsOU 

3 



where /c^ ~ Equation (21) shows that kc is a function of r^. In reality, the 

diffusion leads to the maximum of kc- Therefore tb makes the derivation of the right side of 
Equation (21) be zero. Noting that Vr,b — Vad is approximately proportional to tb — rc, one can 
easily work out the location of point B: 

3 Jc-rs, 



ACsOO ' I 



(22) 



It is found in Fig. 1 that k ^ k Local in the deep convection zone because the turbulent diffusion 
can be ignored here, and the turbulent diffusion dominates in the layer beyond point B. 



-11- 



Using above results, we obtain: 



^ e pcp e V 3 pcp 

Generally, is very small. According to Eq.(20), the length of the adiabatic overshooting 
region is: 



lad ~ e J (24) 

— \- 1 

In the area \r — rc\ < — rd in the overshooting region, the temperature gradient V is 
almost equal to the adiabatic one. In the area \r — rc\ > \rA — rc\, however, according to Eq.(14), 
Eq.(17), and Eq.(5), the temperature gradient V is gradually close to Vr: 



V - V« = (V.. - Vn,A) ■ ^^Pi-\l ^^^tclT^ ^W^^^ ^^^^ 

Although cu in Eq.(24) and Eq.(25) is still unknown, we can estimate it roughly. Equation 
(24) and (25) describe the turbulent motion near the convective boundary, thus we can use co !=:: uc 
where uc is u at the convective boundary. In the deep convection zone, u is almost equal to the 
equilibrium value lJcz = ^ + | which is derived from the localized TCM (see Appendix A). 
u!c < oJcz because the buoyancy is zero at the boundary, and coc > ^ because the diffusion of co. 
Therefore the typical value of uc can be taken as the average, i.e. ujc ~ |(wcz + |)- If Eq.(25) is 



used in the region far away from the convective boundary(beyond the peak of T'T'), oo ^ cue is 
not appropriate. One can use cu = cuo, where Uo is the equilibrium value of u in the overshooting 
region which is introduced in the next subsection. 



Another turbulent correlation is T'T', which can be worked out by using Eq.(l 1): 



T'T' ^0,(\r-rc\ < \rA-rc\) (26) 

And: 

2T T 



T'T' = -TT^i^-d - V)cuk, (|r - rc\ > {ta - rc\) (27) 
-np pg 
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Equation (26) seems to against Cauchy's theorem u[u[T'T' > u'^T' . Actually, T'T' ^ is 



only an approximate solution on the order of {Pe^u'^T'), because Eq.(8) is an approximation on 
that order. Numerical calculations show no confliction. 

Results obtained above are based on C^i = 0. Numerical results of V with both Cei = and 



Cei 7^ are shown in Fig. 2. It is found that the effects of the diffusion of T'T' are only making 



V be smoother. However, there is no adiabatic overshooting region when the diffusion of T'T' 



is present, because T'T' increases near the convective boundary due to the turbulent diffusion 
thus V decreases according to Eq.(8). Numerical results of the turbulent correlations in both 
Cei = and Cei 7^ with different TCM parameters and for different stellar models are shown 
in Figs. 3-5. It is found that the theoretical solutions well fit the numerical solutions in the case of 
Cei = 0. This also validates that the boundary value kc derived from the maximum of diffusion 



is a good approximation. The diffusion of T'T' modifies and smoothes the profile of T'T' and 



u'j.T'. However, k is insensitive to the diffusion of T'T' because that k is mainly dominated by the 



diffusion of itself. The diffusion of T'T' doesn't significantly change the integral value of T'T'. 



According to Eq.(8), the integral value of V or u'^T' is also insensitive to the diffusion of T'T', 
which is found in Figs. (2-5). 



The distribution of T'T' reveals an important property of the overshooting. In the 
nonadiabatic overshooting region, using V ^ Vr, one finds t hat T'T' oc TC^^nd — Vplfc 



according to Eq.(27). This result indicates a maximum of T'T' viXiong 



1985 



Zhang & Li 



2009 ) 



which is shown in Figs. 3-5. Beyond the location of the maximum of T'T', the temperature of a 
turbulent element is gradually close to the temperature of the environment, and the efficiency of 
heat transport significantly decreases. Therefore the area between the convective boundary and 



the location of the maximum of T'T' can be thought as the overshooting region of u'^T'. It is 



found in Figs. 3-5 that the width of the valley of u'^T' is approximately equal to the distance from 



the convective boundary to the location of the maximum of T'T'. In order to get the overshooting 
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length of heat transport, we need to locate the maximum of T'T'. 



Using Eq.(17), defining = S = the decaying index of k (in the case 

of Cei = 0), we get: 

TT ^T{Vad-yR)P'° (28) 



The derivative of T'T' is zero at the peak of T'T'. We get Vr there(denoted as V;^): 

(V^e + ^o)(V,, - V^^) - xV^^ « (29) 

is determined by only one turbulent parameter ^o- 



The typical overshooting length of u'^T' (or V) can be estimated with V\: 
' dlnP ' ~ ' InPc - InP* ' ' InPc - InP* ' 4^ 

Hp 

where Vr^c is Vr at the convective boundary, /v is the distance from the convective boundary to 



the location of the maximum of T'T' and also the typical overshooting length of V. 
/v is worked out as: 

h ~ ^In^Hp (31) 

IXI 

Usually, 1^0 1 is much larger than \x\ and Vad, and V|j can be approximately solved from 
Eq.(29): 

~ (1 - f)Va4 (32) 



Finally, we find: 



= Hk 



(33) 
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where Hk is the scale height of turbulent kinetic energy k defined hy Hk — \ \ . The result 
indicates that V is remarkably modified by the overshooting only in about IHk. It is found in 
Fig. 3 that /y = In^Hk ~ O.Siffc, which is in agreement with Eq.(33). It is shown in Fig.2 that V 
is remarkably modified only in Iff j^. 

3.2. Asymptotic analysis 

In above subsection, we have discussed the turbulent heat transport and the solution of 
turbulent correlations in the overshooting region near the convective boundary based on the 
assumption Cgi — 0. The diffusion of T'T' only modifies turbulent correlations to be smoother 
near the convective boundary. However, it makes more effects on turbulent motions in the 
overshooting region further than IHk away from the convective boundary. In this subsection, we 
investigate the turbulence properties in the outer overshooting region(beyond IHk). 

In the numerical calculations of the TCM, we found that the anisotropic degree u always 
showed an equilibrium value in the overshooting region. A typical numerical result is shown 
in Fig.6. In order to understand it, we discuss the behave of the anisotropic degree u in both 
convection zone and overshooting region, u should be larger than | in the convection zone 
because the buoyancy boosts radial movement of turbulent elements. Actually, to is almost equal 
to the equilibrium value in the convection zone ujcz = 3^ + |(see Appendix A) due to the 
equilibrium between the buoyancy and the return to isotropy term. When turbulent elements go 
across the convective boundary into the overshooting region, the buoyancy prevents convective 
elements moving, thus u decreases to less than | near the convective boundary. However, as u'^T' 
exponentially decreasing, the equilibrium of u) is established again in the overshooting region. 
This results in an asymptotic property of the overshooting region: there is an equilibrium value of 
u! in the overshooting region, ui cjo- 
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By using the asymptotic property a; fa and Approximations 1, 11 & III, it is easy to get the 
asymptotic solution of TCM in the overshooting region(see Appendix B): 

(Ck - l)(uo - \) T kl 

u'rT' = ^ \ ^ -r (34) 



T- 



2 



k^koi^f (36) 

-TO 



where 9 is the asymptotic solution of . 



. = ^1 /^l:- '^- (37) 

a y SC^Wo (1 - Wo) 

k takes the decreasing expression in the overshooting region, which means: '+' is adopted in the 
upward overshooting region and ' in the downward one. 

The equilibrium value ujo is determined by: 

(2C,Ce - C,,CkW + [^C-eilCfc + 2) - C,(Cfc + 2Ce - l)]uJo + \Cs{Ck - 1) = (38) 

The equilibrium value Uo is only a function of turbulent parameters (Cg, Cgi, C^, C^). The 
fact that the buoyancy prevents the radial movement of turbulent elements in the overshooting 
region restricts the turbulent parameters to ensure Uo < \. 

An important thing is where u reaches its equilibrium value Ug. According to Eq.(9), the 
equilibrium of u can be realized only if the buoyancy term synchronically decreases with k 



decreasing. Therefore uj starts to reach its equilibrium value uog beyond the peak of T'T' due to 



\u'^T'\ being decreasing. 

Setting Cgi = in Eq.(38), we find that the asymptotic solution is the same as the results 
in the overshooting region with \r — rc\ > jr^ — rc| by setting u — Uoin Eq(14),(17) & (27). 
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Because Eq.(8) is correct whether Cei = or not, the conclusion that the maximum of T'T' is 
located at about IHk is also correct in both cases. 

It must be mentioned that we have used Approximation I(i.e. Fg ^ 1)» which means that the 
turbulent dissipation is much larger than the thermal dissipation. If k decreases enough to satisfy 
Pe ^ 1, the thermal dissipation should become significant thus T'T' and the turbulent kinetic 
energy k should rapidly decrease to zero. Then uj also rapidly decreases as shown in Fig. 6. In 
another word, turbulent movement can hardly overshoot into the thermal dissipation zone where 

According to discussions above, we can separate the overshooting region into three parts as 
shown in Fig.7: the overshooting region of u'^T' or V with the length of about IH^, the turbulent 
dissipation region in which the asymptotic solution holds, and the thermal dissipation region in 
which the turbulent movement quickly vanishes. The boundaries among those parts are the peak 
of TT and the location ofP^ = l. 



4. Conclusions and discussions 

Turbulent convection models are better tools in dealing with the convective overshooting 
than non-local mixing length theories. However, they are often too complex to be applied in 
the calculations of stellar structure and evolution. In order to investigate the property of the 
convective overshooting and to make it easy to applv turbulen t convection models, we have 



analyzed the TCM developed by Li & Yang (ILi & Yang 



2007h and obtained approximate and 



asymptotic solutions of the TCM in the overshooting region with Pe ^ 1. The main conclusions 
and corresponding discussions are listed as follows: 

1. The overshooting region can be partitioned into three parts: a thin turbulent heat 
flux overshooting region, a power law dissipation region of turbulent kinetic energy, and a 
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thermal dissipation area with rapidly decreasing k. The turbulent fluctuations k, u'^T', and T'T' 
exponentially decrease in the overshooting region as Eqs. (34-36). The equilibrium value of the 
anisotropic degree Uo and the exponential indices of the turbulent fluctuations are only determined 
by the parameters of the TCM. The decaying behaviors of the turbulent fluctuations are similar to 



Xiong & Deng's results (|Xiong 



1989; 



Xiong & Deng 



200lh 



2. The peak of T'T' in the overshooting region is located at about IHk away from the 
convective boundary. In this distance, the modification of V caused by the overshooting is 
remarkable. An approximate profile of V comprises an adiabatic overshooting region with the 
length of lad and an exponentially decreasing function, as described in Eq.(24) and (25). Beyond 
IHk, the modification of V is ignorable and V ~ V/j. It should be noted that the result of IHk 
overshooting distance of turbulent heat transfer is independent of the parameters of TCM, so it 
may be a general property of the overshooting. Our result is similar to Marik & Petrovay(2002) 



whose result shows that the length between the peak of T'T' and the convective boundary is about 
1.2Hk. Meakin & Arnett(2010) simulated the turbulent convection of a 23 Mq star, the data of the 
turbulent kinetic energy and the convective flux in the overshooting region being shown in Fig. 8. 



It is found that the overshooting length of the convective flux u'j.T' is about 0.5 ~ 2Hk which is in 
agreement with our result. 

3. The value of the turbulent kinetic energy at the convective boundary kc can be estimated 
by a method called the maximum of diffusion. The value of turbulent fluctuations in the 
overshooting region can be estimated by using the exponentially decreasing functions and the 
initial value kc- This may significantly simplify the application of the TCM in calculations of the 
stellar structure and evolution. 

There is a distinction between the non-local model of Zahn(1991) and our results, i.e. the 
temperature gradient jumps from nearly adiabatic to radiative in Zahn's model but continuously 
changes in our results (see Fig. 2). This is caused by the assumption in Zahn's model that the 
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turbulent velocity and temperature fluctuation are strongly correlated (|Xiong & Deng 



20011) . In 



u'T' 



our results, the correlativity of turbulent velocity and temperature fluctuation Ryr — . — 

V 2ukT'T' 

quickly decreases to zero then turns to be negative near the convective boundary(see Fig. 9), 
and the asymptotic solution shows that Rvt oc \/~k and exponentially decreases in the turbulent 
dissipation overshooting region. Our result is in agreement with three-dimension simulations such 
as Fig. 6 in Singh et al.(1995) and Fig. 15 in Meakin & Arnett(2007). 
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A. The localized TCM in convection zone. 



The localized TCM results from Eqs.(l-4) by ignoring the d iffusion terms. It is a good 



20011). We attempt to work out the 



approximate of the TCM in the convection zone dLi & Yang 
solution in this appendix. 

Some symbols are defined for conveniences: U = u[T', V = TV, W = Vk, 



Vb-V„ 



Ignoring the diffusion terms of Eqs.(l-4), we get the localized TCM: 

= -—-2BU + 2C,{u--) — 

0= — -BU 



= -2ufAW^ -BV + Ct{l + P, 



(Al) 
(A2) 
(A3) 



, wv 

= -2fAU + 2Ce(l + P; )—^ (A4) 
U = AD{1 - /) (A5) 

Equation (Al) and (A2) show: 

This is the equilibrium value ujcz in convection zone. 
Describing W, V,Uhy f and Pe(= we find: 

CtCeP-\i + p-^y + 2c,uj{i + p-^) + i ^ ^ 

W, V can be worked out as: 

= ABDI{1 - /) (A8) 

According to Pe = Eq.(A8) and Eq.(A7), we get the equation of Pe'. 

aP^ + {b+ 1)P^ + 26P2 + _ at)P, -t^O (AlO) 

where a = 1 + ^ ^ * ^ / determined by / = 1 - ^ according to Eq.(A8). 

Solving Eq.(AlO), we can obtain all turbulent fluctuations of the localized TCM by using 
Eq.(A5), (AS), (A9) and (All). 

An important case is i » 1, thus Pg 1 according to Eq.(AlO). In that case, Eq.(A7) shows: 
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which corresponds to the adiabatic convection. 

Finally, we obtain the turbulent fluctuations according to Eq.(A8), (A5) & (A9): 

fti ABDl (A12) 

U K AD (A 14) 

and the correlativity of turbulent velocity and temperature Rvt'- 



Rvt = , ~ (A15) 



B. Details of deriving the asymptotic solution of the TCM in overshooting region. 

There are the details of obtaining the asymptotic solution of the TCM in overshooting region. 

Some symbols are defined for conveniences: 

A = -^{y - Vad) ~ - - Vad) {Approximation III is used), B = 

Applying the asymptotic property co — cOg and Approximations 1, 11 & III, one can rewrite 
TCM as: 

- (Cfe - i)K - ~ ^^^^ ~ """^ ^^^^ 

d dW^ 

lCsUoj-{W^) = — -BU (B2) 
or or I 

= -BV + 2AujoW'^ (B3) 

d dV C 
ICeiUo-^iW—) ^AU + -^WV (B4) 
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Equation (Bl) is equivalent to: 



(i-.o) -m ^^^^ 



Taking it into Eq.(B2), one gets the equation of W: 

d'W^ 3 , (C,-l)K-i). „ 



ACsUol^ ' {l-ooo) 
Equation (B3) is equivalent to: 



[1 - . \ '' ]W' (B6) 



B 



According to Eq.(B4), (B5) and (B7), one gets another equation of W: 

+ 2CeUJo]W^ (B8) 



d^W^ _ 3 (Cu-l){ujo-\) 



Comparing Eq.(B6) with Eq.(B8), one finds: 
Therefore the equation of ojo is: 

(2C,Ce - CeiCfe)a;,2 + [^Cei(Cfe + 2) - C,(Cfe + 2Ce - l)]ujo + ^C.(Cfc - 1) = (BIO) 



The asymptotic solution of W is derived from Eq.(B6): 



W = W,e.p{±l-, U-ll - '^^ 7/""° - (BID 
2q! V 3G5a;o (1 - ^^o) -H) 



takes the decreasing expression in the overshooting region: '+' is adopted in the upward 
overshooting region and '— ' in the downward one. 
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Fig. 1. — Numerical results of W = \fk, and W Local ~ \j "^^'^^pcp'^"'"'' which is the solution 
of localized TCM (See Appendix A), for the solar model at present age. TCM parameters are: 
a = 0.84, Ck = 2.5, C, = 0.1, Cgi = 0, Cg = 0.2, = 7.0, and Cn = 0.01, but C*, and Cei 
are insensitive to the results. Point C is the boundary of the convection zone, the location of point 
B is calculated by using Eq.(22). 
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Fig. 2. — Numerical results of temperature gradient near the convective boundary in both Cei — 
and Cei 7^ 0, Vo being the temperature gradient of the model with Cei = 0, and Vi corresponding 
to Cei = 0.01. Dotted line Vt, which is almost identical to Vq, is theoretical solution of the 
temperature gradient with Cei — 0. The stellar model and other TCM parameters are the same as 
Fig.l. Point A is the boundary of the adiabatic overshooting region. Our theoretical result shows 
lad ~ O.OlSHp in those TCM parameters set, the numerical calculation being 0.015Hp. 
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Fig. 3. — Numerical results of T'T', u'^T', k near the convective boundary in both Cgi = and 
Cei ^ 0, where U = i^T, W = Vk,V = TT. Dashed lines correspond to C^i = 0, soUd Unes 
to Cei = 0.01. Dotted lines are the theoretical solutions with Cgi — 0. The stellar model and other 
TCM parameters are the same as Fig.l. 
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Fig. 4. — Numerical results of T'T' , u'^T' , k near the convective boundary in both Cgi = and 
Cei 7^ 0, where U = ^T', W ^ \/k,V ^ TT. Dashed lines correspond to Cei = 0, soUd lines 
to Cei = 0.01. Dotted lines are the theoretical solutions with Cei = 0. The stellar model is a 7Mq 

star model at the top of RGB phase. Others TCM parameters are: a = 1.0, C^ = 2.2, C^ = 0.1, 
Ce = 1.0, and Q = 4.0, Ca = 0.01. 
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Fig. 5. — Numerical results of T'T', u'^T', k near the boundary of the convective core in both 
Cei = and C^i ^ 0, where U = ^F, W ^ Vk,V ^ TT. Dashed lines correspond to 
Cei = 0, solid lines to C^i — 0.01. Dotted lines are the theoretical solutions with Cgi = 0. 
The stellar model is an early main sequence model of a SMq star. Others TCM parameters are: 
a = 1.0, Ck = 2.1, Cs = 0.2, Ce = 0.5, and Ct = 3.0, Cn = 0.01. 
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Fig. 6. — Numerical result of the structure of co in overshooting region. The stellar model is the 
solar model at present age. Cgi = 0.01. The others TCM parameters are the same as Fig.l, except 
a = 0.2 in order to enlarge 6 to show the thermal dissipation region in which Pg <S 1. With those 
parameters, the equilibrium value in convection zone is Ucz = 0.6, and the equilibrium value in 
overshooting region is cOg = 0.293 which denoted as the dotted line. 
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Fig. 7. — The structure of the overshooting region. K = k, U = u[.T', V = T'T'. The stellar 
model is the solar model at present age. C^i = 0.01, the others TCM parameters are the same as 
Fig.l, except a — 0.2. With those parameters, in the turbulent dissipation region with P^. » 1, 
theoretical result shows 9 — 17.5 vs the numerical result 17.6, theoretical result of exponential 
decreasing index of U"^ being 26.3 vs the numerical result about 25.6. K is almost parallel to V, 
which is in consistent with the asymptotic solution. In the thermal dissipation region with Pe "C 1> 
turbulent motion vanishes. 
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Fig. 8. — Numerical data of Meakin & Amett (2010)'s results. The data of model 'hi' in their 

paper are plotted, where U = Fc = pCpu[.T' . Only the downward overshooting region is shown. 
The distance from the convective boundary (where u'^T' — 0, about R — 0.62 x lO^cm) to the 
right part of the valley of u'^T' is about 0.5 ~ 2iffc. 
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Fig. 9. — Numerical results of the correlativity of turbulent velocity and temperature Ryr- The 
stellar model is the solar model at present age. Other TCM parameters are the same as Fig.l. Ryr 
rapidly decreases to zero in the overshooting region. In the convection zone near the convective 
boundary, the diffusion significantly enlarges T'T' when Cei ^ (see Fig.3), and then Rvt is very 
small. In the interior of convection zone, localized TCM shows the equilibrium value of Ryr is 
RvT,cz = \l '^%Z^'c^ (^^^ Appendix A). The TCM parameters show Rvt,cz = 0.384. 



